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The temperature dependence of spin currents in insulators at the finite temperature 
near zero Kelvin is theoretically studied. The spin currents are carried by Jordan- Wigner 
' (^ ' fermions and magnons in one- and three- dimensional insulators. The quasiparticle de- 

^Z^ ' scription of one-dimensional spin systems is valid only in the finite temperature near 

O , zero Kelvin. These spin currents are generated by the external magnetic field gradient 

O 1 along the quantization axis and also by the two-particle interaction gradient. In one- 

dimensional insulators, quantum fiuctuations are strong and the spin current carried by 
Jordan- Wigner fermions shows the stronger dependence on temperatures than the one 
1...^ ^ by magnons. 

i^ ! 
l 1 Keywords: Spin currents; Jordan-Wigner fermions; Magnons; Insulators; Spintronics. 

OO 

*^ \ 1. Introduction 

\mm^ , Recently a new branch of physics and nanotechnology cahed spintronica^J^ has 

emerged and spin transport phenomena in condensed matter systems have been 
attracting special interests because of applications to spintronics. This research field 
has seen a rapid development over the last decades.!^ The aim of spintronics is the 
control and utilization of the spin as well as charge degrees of freedom of electrons. 
Spintronics avoids the dissipation from Joule heating by replacing charge currents 
C^ \ with spin currents, and therefore to clarify the properties of spin currents is an 

important t heor etical issue from viewpoints of fundamental science and potential 
applications v^^ 

The spin current means a flow of the spin angular momentum, in general. It 
also flows in insulators as well as in metals where conduction electrons carry a spin 
current. In insulators there is no conduction electrons, but there exists an other 
kind of carrier, namely, spin- waves. 

Experimentally, a spin-wave spin current, a spin current carried by spin-waves 
has already been established as a physical quantity. Y.Kajiwara et al.'^have reported 
that a spin-wave spin current in an insulator can be generated and detected using 
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direct and inverse spin-Hall effects!^ Moreover through the effects, it is possible to 
convert an electric signal in a metal into a spin- wave spin current, and vice versa. 
This spin-wave spin current has a novel feature that this current persists for much 
greater distance than the conduction electron spin current, which disappears within 
very short distance (see Table [1] ) . This novel feature of the spin- wave spin current 
is expected to lead to developments beyond silicon-based technologies. 

Theoretically, Meier and Loss^ have studied the magnon transport in both fer- 
romagnetic and antiferromagnetic materials. They have found that the spin conduc- 
tance is quantized in the units of order (g/ie)^/^ in the antiferromagnetic isotropic 
materialsfg; the gyromagnetic ratio, /ib; the Bohr magneton, h; Planck constant). 
Fujimoto^ has studied a spin Hall effect of spin- waves in frustrated magnets. A 
longitudinal magnetic field gradient induces a transverse spin current carried by 
spin wave; Jx — a^^^dB/dy, c^™ oc T^ at the low-temperature, where J^ is the 
spin Hall current, cr^?^ is the Hall conductivity for the spin current and B is the 
magnetic field along the quantization axis. Katsura et al! I I ^have reported the 
intrinsic thermal Hall effect for magnons due to the anomalous velocity. 

Though thus theoretical studies of spin currents in insulators have developed, it 
should be emphasized that, besides magnons, there exists an other kind of carrier in 
insulators, namely, the Jordan- Wigner (J-W) fermion. This carrier is peculiar to one- 
dimensional spin systems. In order to clarify the properties of spin current s carried 
by J-W fermions, we focus on the one-dimensional spin-1/2 XXZ model! I I I 
The model reads 

Hxxz = rY^iSfSf^, + SfSf^, + A,5f 5f+i), (1) 

i 

where Sf is a component of a spin-1/2 operator on / th site {a — x,y,z). The 
constant parameter F is the exchange coupling constant, and A is the anisotropic 
interaction parameter. This model can be well mapped into spinless fermion systems 
via the J-W transformation. When A = 1 the XXZ model reduces to the Heisenberg 
model, and when A = this model reduces to the XY model; Hxy = P J2i i'SfSf^^ + 
SfSy_^_^), which describes the disorder phase; {S") — 0. After a transformation which 
respects the spin commutation relations; Sf — >■ (— l)*5f, Sf -^ {—lySf, Sf — )■ S^, 
the XXZ model changes the sign of parameters as F — )• — F, A ^- A. Thus regardless 
of the signs of the parameters, it is enough to consider the case F > 0. This paper 
focuses on the properties of the spin current when the value of | A | is smaller than 
one. 

Quite recently, theoretical studies for spin currents i n o ne-dimensional spin 
chains also have been rapidly progressing. Trauzettel et al\^ have calculated the 
ac magnetization current and the power absorption of the XXZ model. Hoogdalem 
et al\^ have considered rectification effect s ■^'^'-■'^' through bosonization techniques 
in both ferromagnetic and antiferromagnetic systems. 

In this paper, we study the effects of external magnetic field gradients along the 
quantization axis and two-particle interaction gradients on the spin current carried 
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by J-W fermions in one-dimensional insulators; the quasiparticle description ofone- 
dimensional spin systems is valid in the finite temperature near zero Kelvin! I I 
The magnetic field is equivalent to a chemical potential for J-W fermions and the 
gradient plays a role similar to an electric field in electron systems!^ The two- 
particle interaction arises from the anisotropic interaction in eq. ([1]) and the gradient 
also generates the spin current. The main aim of this paper is to reveal the temper- 
ature dependence of the spin current carried by J-W fermions and then compare it 
with that of the magnon current in the finite temperature near zero Kelvin. 

This paper is structured as follows. In the one-dimensional XXZ model, the spin 
degrees of freedom can be mapped efficiently into fermion degrees of freedom via the 
J-W transformation. Thus first, we define microscopically the spin current density 
carried by J-W fermions. Second, we evaluate it through the standard procedure 
of the Schwinger-Keldysh Green's function. The effects of magnetic field gradients 
and two-particle interaction gradients on the spin current are clarified in SeclJJ 
Then the magnon current generated by the two-particle interaction (i.e. magnon- 
magnon interaction) is also calculated in Sec|3l Last, the temperature dependence 
of the spin current carried by Jordan- Wigner fermions is compared with that of 
the magnon current in Sec|4l We discuss the cause for the stronger dependence of 
the spin current carried by J-W fermions on temperatures than that of the magnon 
current in the finite temperature near zero Kelvin. 

Table 1. The features of spin currents; spin currents carried by conduction electrons (i.e. conduction electron 
spin currents), by magnons (i.e. magnon currents) and by J-W fermions (i.e. J-W spin currents). The magnon 
current, which is the quantized spin-wave spin current, persists for much greater distance than the conduction 
electron spin current because the dissipation from Joule heating does not exist in insulators. 

Features Conduction electron spin current Magnon currents (3-dim) J-W spin currents (1-dim) 



• Systems 


Metals 


Insulators 


Insulators 


• Spins 


Conduction electrons' spins 


Localized spins 


Localized spins 


• Statistics 


Fermions 


Bosons 


Fermions 


• Decay length 


A few nano-meters 


A few centi-meters - 


\ 



2. Jordan- Wigner Spin Current 

The Mermin- Wagner theorem states that continuous symmetries cannot be spon- 
taneously broken at the finite temperature in one- and two-dimensional systems. 
Therefore in such low-dimensional spin systems, the Holstein-Primakoff transfor- 
mation is not useful. In other words, because the one-dimensional XXZ model has 
SO (2) symmetry that the magnon description cannot be well applied. Thus in one- 
dimensional spin systems, the spin degrees of freedom cannot be mapped efficiently 
into boson degrees of freedom, but it can be well mapped into spinless fermion de- 
grees of freedom via the J-W transformation. Then the XY model reduces to a free 
fermion system. 
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2.1. Definition 

The J-W transformation maps the spin operators of the XXZ model onto fermionic 
operators; S+ = cjU^Hl ~ 2nj),Sr = c,n}zi(l - 2n,), 5f ^ n, - 1/2, n, = c\c,. 
Operators c^ /c are J-W fermion creation/annihilation operators which satisfy the 
fermionic commutation relation. The operator, n = c^c, is the number operator 
of the J-W fermion. This spinless fermion carries the spin angular momentum of 
localized spins in one-dimensional systems. The z-component (i.e. S^ — n^ ~ 1/2) 
shows that a spin down can be viewed as an empty site and a spin up corresponds 
to the presence of a fermion. After the J-W transformation and the canonical trans- 
formation; a — > (— l)'ci, T-Lxxz reads 

-Hxxz = -^^(c|c,+i +c|+ic,) +r^ A,(n, - i)(n,+i - i). (2) 

i i 

Eq.([2]) shows that only fcrmions between neighbors feel an interaction FA and 
J-W fermions hop between neighboring site with a hopping matrix element F/2. 
Therefore in the continuous limit, Hxxz reads Hq + Va as an effective Hamiltoniaro; 

r 1 (f 

Ho = J dxcHx,t){- — — )c{x,t), (3) 

Va = dxJ{x)c''{x,t)c^x,t)c{x,t)c{x,t). (4) 

Here the parameter m is the effective mass of a J-W fermion which is related to the 
curvature of the dispersion relation, and the parameter J corresponds to —FA in 
the discrete model. 

Moreover, there are magnetic impurities in real materials. Magnetic impurity 
scatterings make the lifetime of quasiparticles r, which is inversely proportional to 
the imaginary part of the self-energy, finite. Here it should be emphasized that the 
lifetime of quasiparticles in insulators is generally far longer than that of conduc- 
tion electrons in metals because the dissipation from Joule heating does not exist in 
insulators; the lifetime of conduction electrons in metals is mainly caused by lattice 
defects, nonmagnetic impurities and magnetic impurities et al. through Coulomb 
interactions. Because quasiparticles in insulators have no charge degrees of freedom 
that the infiuence of impurities is far smaller than that in metals; magnetic impuri- 
ties, at most, might cause impurity scatterings. Furthermore in real materials, the 
rate of impurities such as lattice defects and nonmagnetic impurities is, in general, 
far larger than that of magnetic impurities. Therefore the imaginary part of the 
self-energy due to magnetic impurity scatterings can be assumed to be so small 
that the dispersion relation does not change (i.e. ujk oc fc^), and we omit the vertex 
correction (i.e. the diffusive spin current -). Let us mention that in eq.Q we have 
defined the effective mass of a J-W fermion, m, which includes the contribution of 
the self-energy due to magnetic impurity scatterings. 



^Through this paper, we take h = 1 for convenience. 
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Thus we introduce the hfetime as a phenomenological parameter and assume 
that it is independent of temperatures. In order to take account of magnetic im- 
purity scatterings, we adopt retarded and advanced Green's functions, G\^ = 
[uj — Wk + i/(2T)]~^ = {G\^)* where uJk = fc^/(2m), on perturbative expansion 
in SubsecH^] and Subsec l^31 Through these procedure,!^ the effects of magnetic 
impurity scatterings, which are the essential dynamics in real materials, are phe- 
nomenologically included into our calculation. 

The density of J-W fermions , pjwj of the system is defined as the expectation 
value of the number operator of J-W fermions as 

P}w{x,t) = {c\x,t)c{x,t)). (5) 

Through Heisenberg's equation of motion, the spin current density carried by J-W 
fermions, i.e. the J-W spin current, jjwj is defined as 

dtpjw + dxjjw = 0, (6) 

jjw = -Re[i((9.ct)c)]. (7) 



2.2. Evaluation: magnetic field along the quantization axis 

In this subsection, we consider the J = case and evaluate the J-W spin current 
generated by magnetic field gradients along the quantization axis. We treat Vb! 
Vb — ~ J dxB{x)S^ — ~ J dxB{x)[c' {x,t)c{x,t) — 1/2], as a perturbative term to 
evaluate the first-order contribution in B without vertex corrections. 

Through the standard procedure of the Sch winger-Ke ldysh (or contour-ordered) 
Green's function, B5 l ^6 | 27 | ^j^g Langreth method,'23I2SE91 the J-W spin current jf^ 
is evaluated as 

The variable L is the chain length and the variables G*, G^, and G^ are the time- 
ordered, greater, and lesser Green's functions, respectively. Fermionic greater and 
lesser Green's functions are defined as G>(ii,i2) = —i{c{ti)c'{t2)),G'^{ti,t2) = 
i{cHt2)c{t^)). They satisfy the relations: G*^„ = Gl^ + G>^, Gl."G<. = G^^- 
Gt u] ■ Thus the current is rewritten as 



J% = -2^^I°^E.,(fc + §)e-«^S- Jdo; 






'\GkJ^q/2,LoGk-q/2,u} + ""fe+g/2,w""fe~g/2,w + '-'A:+q/2,w'-'/£-ij/2,w )■ \^) 

The retarded Green's function, G'', is given by G\^^ = [uj-ujk + i/{2T)]^^ = {G%^^)* 
and T is the lifetime of a J-W fermion. The energy a;^ is oj^ = Fk^, F = {2m)^^. 
Then the function G' reads G'j^_^/2,^ = Gl^^ + {GlJ'^{-Fkq + Fq^/A) + 0{{kqf). 
This approximation demands that the magnetic field varies slowly and moderately 
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in space (compared to the J-W fermion mean- free path), and the value of dxB is a 
constant. Then eq.® reads 

•| [4t(1 - Hk) + i/3e'5(-^'='^'^)nfc2]Ffcg 

-[4t(1 - uk) - i^e^f^'^'-^^nfc^jFfcg}, (10) 

where n^ is the Fermi distribution function, fi is the chemical potential, /3 = l/{kBT) 
and fee is the Boltzmann constant. We have used the approximation; nj^^g/2 ~ 
Uk - pPe^'^^^^^i-'-'^nlkq. By the Sommerfeld expansion (T < T^.T^ ; the Fermi 
temperatur^j) , the chemical potential is evaluated as /i(r) = ep + (7''fcB7')^/(6eF) + 
0{T'^), where ep is the Fermi energy. 

Finally, the J-W spin current generated by the magnetic field gradient along the 
quantization axis in the finite temperature near zero Kelvin reads 



R 4/Fr,, vr^ / (7rfcBr)2 ^ , , 



2.3. Evaluation: two-particle interaction 

In this subsection, we treat the B = Q case and evaluate the J-W 
spin current generated by the two-particle interaction gradient. The two- 
particle interaction arises from the anisotropic exchange interaction Va; Va = 
j dxJ{x)c^ {x^t)c\x,t)c{x^t)c{x^t), and we treat it as a perturbative term to eval- 
uate the first-order contribution in J without vertex corrections. 

Through the same procedure with the last subsection, the J-W spin current 
generated by the two-particle interaction Va is evaluated as 

'\^k+nl2.u^k-al2.uj ^ ^k+al2.i^^ k-al2M ^ ^k+al2.L0^k 



^k+q/2,i^^k-q/2,u} ^ ^ k+q/2,LU^ k-q/2,uj ^ ^ k+q/2,i^^ k-q/2,uj 

+ ^k+q/2,u:^k-q/2,iJ ^ ^k+q/2,ij^k-q/2,uj) + ^i-J )• (12) 

We have used the relation that the numbe r operator, {c^{t)c{t)), is the same as 
the equal-time lesser Green's function^^l^ {c''{t)c(t)) = —iG'^{t,t). Because each 
Green's function is not independent; G^ — G^ — G^ — G^, that eq. ([T2l) is rewritten 
as 



l^^L^ '-'-'-'- L^kqy"- ' 2'^ "^-l" L^p J ""^ ^p,uj' 

'y-'k+q/2,uGk-q/2,uj + ^ k+q/2,uG k-q/2,uj + ^fc+<j/2,w^fe-g/2,w)- \^'^) 

Also in this case, we assume that the two-particle interaction varies slowly and 
moderately in space (compared to the J-W fermion mean- free path), and the value 



The Fermi temperature of this system(insulator) is estimated as a few hundred Kelvin. 
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of dxJ is a constant. Therefore the retarded green's function can be approximated; 
Gl-^/2,u. = Gl^ + (GIJH-Fkq + Fq^A) + 0{{kqf). Then eq.dni) reads 

3% - -i^Re5:/3fe^(^'=^--)fc^n^^n,a.J. (14) 

k p 

Finally, the J-W spin current generated by the two-particle interaction in the finite 
temperature near zero Kelvin reads 



4t „ „ TT 



2 



[Tlk^Tfl {Tlk^Tf 



Jtw = -^(1^2 + ^)<! 2 6p + ^^^^ r ""Tl.^.1 >9.J- (15) 



Gep 



12 



^F + ^^HP^ 



3. Magnon Current 

According to the Mermin- Wagner theorem, three-dimensional spin systems are the 
most suitable set-ups for the study of magnon transports. Therefore in order to 
clarify the thermal properties of magnon currents, we focus on three dimensional 
ferromagnetic insulators at the low-temperature near zero Kelvin, which has no 
conduction electrons. The spin degree of freedom reduces to the magnon degree of 
freedom via the Holstein-Primakoff transformation; the transverse components of 
the exchange interactions give rise to hopping of the magnons, while the longitudinal 
component give rise to the interaction. 

External magnetic fields, of course, generate magnon currents also in this case; 
the magnetic field along the quantization axis is equivalent to a chemical potential 
for magnons as well as J-W fermions, and the gradient generates magnon cur- 
rents. Moreover time-dependent transverse magnetic fields cause quantum fiuctua- 
tions and generate a magnon current through pumping effects.'^ Here, however, we 
consider three dimensional ferromagnetic insulators without any external magnetic 
fields. 



3.1. Definition 

The low-energy effective HamiltoniarF^ for magnons in three-dimensional insulators 
reads Hm = 'Ho(m) + Kn; 

Ho(m)= /"d3a;at(x,t)(-^-^)a(x,i), (16) 



Vm= / d^a:;J(m)(x)at(x,t)a1'(x,i)a(x,i)a(x,i). (17) 

Here mmag is the effective mass of a magnon, and 1^ represents the magnon-magnon 
interaction. Operators a^ /a are magnon creation/annihilation operators which sat- 
isfy the bosonic commutation relatio n. In the low-temperature, only the two-particle 
interaction J(iji) is important! II 
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The magnetic impurity scattering is, of course, an inevitable dynamics in real 
materials even when the system is under low-temperature. It makes the lifetime 
of magnons, Tm, finite. Therefore by the same procedure with Subsec l2Tl we phe- 
nomenologically include the effects of magnetic impurity scatterings, which are the 
inevitable dynamics in real materials, into our calculation. Also in this case, we 
have assumed that the imaginary part of the self-energy due to magnetic impurity 
scatterings is so small that the dispersion relation does not change, and we omit 
the vertex correction (i.e. the diffusive spin current). Furthermore, we assume that 
the lifetime of magnons is independent of temperatures. 

The magnon density, pm, of the system is defined as the expectation value of 
the number operator of magnons 

p„,(x,i) = (at(x,i)a(x,i)). (18) 

Through Heisenberg's equation of motion, the magnon current (density), jm, is 
defined as 

dpm 
dt 

Then the magnon current reads 



V-j,n = 0. (19) 



j;;(x,t) = ^Re[i((a,at(x,t))a(x,i))], (20) 

where i/ is a direction for a magnon current to flow (v — x,y, z). 

3.2. Evaluation 

From now on, we treat the magnon-magnon interaction term, T4n, as a perturba- 
tive one to evaluate the first-order contribution in Jj^) without vertex corrections. 
Through the same procedure with the last section, the Langreth methodj I I I 
the magnon current is evaluated as 

'('-'k+q/2,w^k-q/2,w + '-'k+q/2,w'-'k-q/2,w + '^k+q/2,w^k-q/2,Lj) + ^("^(m) X^l) 

Here ^ is a volume of the system. Bosonic greater and lesser Green's functions 
are defined as G^{ti,t2) = —i{a{tx)a^{t2)),G^{ti,t2) = —i{a^{t2)a{ti))] the sign of 
the lesser green's function is different from the fermionic case. The retarded Green's 
function is GJ^ j^ = [w— Wk+*/(2T„i)]^^ = (G^ ^)*, where the energy (Dk is Wk = -Dfc^, 
D = l/(27Tiniag) and Tni is the lifetime of magnons. The lifetime represents the 
damping of spins (rm is inversely proportional to the Gilbert damping parameter, 
a). Then the retarded Green's function satisfies Gj^+q/2 cj — ^k w + (^kt>j)^(-^^ ' 
q + q^/4) -I- 0((k • q)^). This approximation demands that the magnon-magnon 
interaction varies slowly and moderately (compared to the magnon mean- free path), 
and the value of dvJ(m) is a constant. 
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Thus eg. d^Tt is rewritten as 

T77 



^^^ReEk,q(^. + f )e-'''-^JM-q ^ C(i)r(i)(¥)'/' 



•{[/3(e"^^'=' + 26^2^^'=') + 4ir„,(l + nk)]k • q 
+ [/3(e"^^'=' + 20-2/5^'^') - AiT^[l + nk)]k • q} (22) 

Here fik is the Bose distribution function, Q is the Riemann zeta function, C(3/2) = 
2.612, and V is the Euler gamma function, r(3/2) = \/tt/2. Because we consider the 
the low-temperature regime near zero Kelvin, we have approximated as rik+q/2 ~ 
nk - l3Dc-^^^^ (1 - e-'3^'=')-2k • q « % - pD{e-^^^^ + 2e-2/5£'fc')k . q. Moreover 
we have used the isotropy condition in the direction for k; ^j. fc^k • q = ^j^. fc^Qj^/S. 
Finahy in the low-temperature near zero Kelvin, the magnon current generated 
by the magnon-magnon interaction gradient reads 



J in 



(i + ^)C(^)r(^)(^)^5.^M. (24) 



Table 2. The temperature dependence of J-W spin currents and magnon ones generated by the 
magnetic field gradient and by the two-particle interaction gradient. The constant 7 is defined as 
7 = 6[ep/{7rfcB)]^. J-W spin currents exhibit the stronger dependence on temperatures than magnon 
ones in the finite temperature near zero Kelvin (i.e. ^ T < 1 [K]). 

Temperature dependence J-W spin currents (1-dim) Magnon currents (3-dim) 

• Magnetic field; d^B oc ^fpPf^ oc T^/^ (Refl6ll34)l 

• Two-particle interaction; dxJ, (dvJ(m)) oc [7 -|- T^ - t2/{4 + iT"^ h)] oc T^ 



4. Temperature Dependence near Zero Kelvin 

The Mermin- Wagner theorem gives crucial differences between low-dimensional spin 
systems and three-dimensional ones. It prohibits the spontaneous breaking of the 
continuous symmetry of the low-dimensional XXZ model, i.e. S0(2) symmetry, 
in the finite temperature. In other words, in low-dimensional systems quantum 
fluctuations are so strong as to destroy magnetic orders. Then the thermal properties 
of each spin current in insulators, the J-W spin current and the magnon one, are 
drastically different from each other in the finite temperature near zero Kelvin 
(i.e. ^ T < 1 [K]); as shown in Table [21 the J-W spin current exhibits the 
stronger dependence on temperatures than the magnon current, i.e. the exponents of 
temperatures in J-W spin currents are smaller than those in magnon currents. Thus, 
the influence of quantum and thermal fluctuations in one-dimensional insulators is 
larger than that in three-dimensional ones. The Table [2] also shows that the current 
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generated by the two-particle interaction gradient shows the weaker dependence on 
temperatures than the one by the magnetic field gradient along the quantization 
axis. 

5. Summary and Discussion 

We have theoretically studied the temperature dependence of the spin current in 
one- and three-dimensional insulators. In one-dimensional insulators, the spin cur- 
rent is carried by Jordan- Wigner fermions. In this system, quantum fluctuations are 
strong and the J-W spin current shows the stronger dependence on temperatures 
than the magnon one in the finite temperature near zero Kelvin. 

Experimentally though the direct measurement of spin currents is impossible 
at this stage, the spin currents we have discussed would be identified by observing 
the temperature dependence of the spin currents through inverse spin-Hall effects. 
Theoretically, the microscopic calculation of the lifetime caused by magnetic impu- 
rity scatterings and also by phonon ones et al., and the diffusive spin current with 
vertex corrections is an significant theoretical issue. 

Acknowledgements 

The author would like to thank K.Totsuka, G.Tatara, S.Fujimoto and M.Oshikawa. 
He is also grateful to Y.Korai and K.Taguchi for invaluable discussion and T.Kimura 
for reading of the manuscript and useful comments. We wish to acknowledge the 
G.Tatara-laboratory members in Tokyo Metropolitan University for warm hospital- 
ity during his stay. 

This work was supported by the Grant-in-Aid for the Global COE Program 
"The Next Generation of Physics, Spun from Universality and Emergence" from 
the Ministry of Education, Culture, Sports, Science and Technology (MEXT) of 
Japan. 

References 

1. S.Maekawa, Concepts in Spin Electronics (Oxford science publications, 2006). 

2. Y.Tserkovnyak, A.Brataas, G. E.W.Bauer, and B.I.Halperin, Rev.Mod.Phys. 77, 1375 
(2005). 

3. D.D.Awschalom and M.E.Flatte, Nature Phys 3, 153 (2007). 

4. S. A. Wolf, D. D. Awschalom, R. A. Buhrman, J. M. Daughton, S. von Molnar, M. 
L. Roukes, A. Y. Chtchelkanova, and D. M. Treger, Science 294, 1488 (2001). 

5. I.Zutic, J.Fabian, and S.D.Sarma, Rev.Mod.Phys. 76, 323 (2004). 

6. Y.Kajiwara, K.Harii, S.Takahashi, J.Ohe, K.Uchida, M.Mizuguchi, H.Umezawa, 
H.Kawai, K.Ando, K.Takanashi, S.Maekawa, and E.Saitoh, Nature 464, 262 (2010). 

7. Y.K.Kato, R.C.Myers, A.C.Gossard, and D.D.Awschalom, Science 306, 1910 (2004). 

8. F.Meier and D.Loss, Phys. Rev. Lett. 90, 167204 (2003). 

9. S.Fujimoto, Phys. Rev. Lett. 103, 047203 (2009). 

10. H.Katsura, N.Nagaosa, and P.A.Lee, Phys. Rev. Lett. 104, 066403 (2010). 

11. Y.Onose, T.Ideue, H.Katsura, Y.Shiomi. N.Nagaosa, and Y.Tokura, Science 329, 297 
(2010). 



August 1, 2011 0:8 WSPC/INSTRUCTION FILE IJMPB0708 



Spin Currents in Insulators 11 

12. R.Matsumoto and S.Murakami, Phys. Rev. Lett. 106, 197202 (2011). 

13. E.Fradkin, Field Theories of Condensed Matter Systems (Allan M.Wylde,1991). 

14. I. Affleck, in Fields, Strings and Critical Phenomena, Les Houches Summer School 
Proceedings, Vol 49 ( Elsevier Science, Amsterdam, 1989). 

15. T.Giamarchi, Quantum Physics in One Dimension ( Oxford University Press, 2003). 

16. B.Trauzettel, P.Simon, and D.Loss, Phys. Rev. Lett. 101, 017202 (2008). 

17. K.A.van. Hoogdalem and D.Loss, "Rectification of spin currents in spin chains", 
'cond-mat/1102480L. 

18. D.E.Feldman, S.Scheidl, and V.M.Vinokur, Phys. Rev. Lett. 94, 186809 (2005). 

19. B.Braunecker, D.E.Feldman, and F.Li, Phys. Rev. B 76, 085119 (2007). 

20. A.Virosztek and J.Ruvalds, Phys. Rev. B 42, 4064 (1990). 

21. S.Fujimoto, Bussei Kenkyu. 54, 207 (1990). 

22. F.D.M.Haldane and D.P.Arovas, Phys. Rev. B 52, 4223 (1995). 

23. A.Takeuchi, K.Hosono, and G.Tatara, Phys. Rev. B 81, 144405 (2010). 

24. H. Haug and A. P. Jauho, Quantum Kinetics in Transport and Optics of Semiconduc- 
tors (Springer New York, 2007). 

25. J. Rammer and H.Smith, Rev.Mod.Phys. 58, 323 (1986). 

26. A.Kamenev, "Many-body theory of non-equilibrium systems", cond-mat/0412296[ 

27. T.Kita, Prog.Theor.Phys. 123, 581 (2010). 

28. G.Tatara, H.Kohno, and J.Shibata, Physics Report. 468, 213 (2008). 

29. D. A. Ryndyk, R. Gutierrez, B. Song and G. Cuniberti, Energy Transfer Dynamics 
in Biomaterial Systems, (Springer- Verlag, 2009), pp. 213-335. 

30. K.Nakata and G.Tatara, J. Phys. Soc. Jpn. 80, 054602 (2011). 

31. T.Nikuni, M.Oshikawa, A.Oosawa, and H.Tanaka, Phys. Rev. Lett. 84, 5868 (2000). 

32. A.K.Kolezhuk, Phys. Rev. B 59, 4181 (1999). 

33. T.Giamarchi and A.M.Tsvelik, Phys. Rev. B 59, 11398 (1999). 

34. K.Taguchi and G.Tatara, unpublished. 



